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a^ : ABSTRACT 
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(N 

O ' A theoretical description is presented for low-temperature magnetic-field induced three- 

dimensional (3D) ordering transitions in strongly anisotropic quantum antiferromag- 

Q\ ', nets, consisting of weakly coupled antiferromagnetic spin-1/2 chains and ladders. First, 

effective continuum field theories are derived for the one-dimensional subsystems. Then 
the Luttinger parameters, which determine the low-temperature susceptibilities of the 
chains and ladders, are calculated from the Bethe ansatz solution for these effective 
models. The 3D ordering transition line is obtained using a random phase approxima- 



O ! tion for the weak inter-chain (inter-ladder) coupling. Finally, considering a Ginzburg 

^ ' criterion, the fluctuation corrections to this approach are shown to be small. The na- 

ture of the 3D ordered phase resembles a Bose condensate of integer-spin magnons. It 
^ ■ is proposed that for systems with higher spin degrees of freedom, e.g. N-leg spin-1/2 

j^ . ladders, multi-component condensates can occur at high magnetic fields. 



1 Introduction 

Compounds of weakly coupled spin chains typically have an ordering transition from 
a high-temperature quasi-one-dimensional (ID) phase to a low-temperature three- 
dimensional (3D) phase at a critical temperature which depends on the inter-chain 
coupling constant |jl|. This transition can be suppressed if the individual chains are 
gapped spin liquids, as it is the case for Haldane or spin-Peierls systems, or for com- 
pounds with an Ising-anisotropy in the intra-chain exchange coupling. In these systems, 
the low-energy spectrum consists of a singlet ground state with an excitation gap to 
the first triplet. This gap can be reduced and eventually overcome by turning on and 
increasing an external magnetic field. Once the spin gap is destroyed, the residual inter- 
chain coupling can lead to 3D ordering at low temperatures. 0, H, Q In this paper, we 
discuss a quantitative theoretical approach to study such magnetic-field induced tran- 
sitions, based on an exact field-theoretical description of the low-energy intra-chain 
dynamics which drives the transition, combined with a mean-field theory (including 
quantum fluctuation corrections) for the inter-chain exchange. 

Our results are in good agreement with recent experiments on the compounds 
TlCuCls^ and CM2(C2-ffi2^2)2C'^4 0, where transition lines hc(T) were extracted 
from an analysis of the temperature-dependent magnetization and from NMR data 
for the 1/Ti relaxation. Typically, the spin gap in most of the ladder compounds 
known to date is too large to be overcome by presently available magnetic fields. How- 
ever, these particular materials have small spin gaps of the order 10 - 20 Kelvin, which 
makes the interesting gapless regime experimentally accessible. It has recently been 
pointed out that Cu2{C2Hi2N2)2CU may better be modeled as an ensemble of weakly 
coupled dimers than as an antiferromagnetic 2- leg ladder 0]. Whatever the precise 
structure may turn out to be, a magnetic-field induced ordering transition can occur in 
all anisotropic spin systems with a singlet-triplet excitation gap, including weakly cou- 
pled Ising-like chains, spin-Peierls chains, and ensembles of spin dimers. Other possible 
candidate materials with spin gaps include KCuCl^ [Bl, CuGeO^ M, a' — NaV205\^^, 



and the homologous series of cuprates SrnCun+i02n+i [0- 

The exact nature of these 3D ordered phases is currently under debate [Q §]. In 
the case of weakly coupled gapless chains which undergo a 3D ordering transition 
even in the absence of a magnetic field, long-range antiferromagnetic order is found 
below the transition temperature |I[]. For the spin-gapped compounds TlCuCl^ and 
Cu2{C2Hi2N2)2CU the field-induced low-temperature transition resembles that of a 
Bose-Einstein condensation of integer-spin magnons. The critical exponent for the 
transition line of such a condensate, i.e. he oc T" with a = 3/2, is rather close to the 
experimentally observed behavior |Q, ^. 

In this work a numerical solution of the Bethe ansatz is used to determine the 
susceptibilities of the one-dimensional subsystems, combined with a generalized RPA 
approach for the low-temperature 3D ordering transition. In the following section we 



illustrate this approach by discussing the case of weakly coupled Heisenberg chains 
with an easy-axis anisotropy (XXZ model). Subsequently, corrections due to fluctua- 
tion effects are determined, which turn out to be rather small. Then the case of weakly 
coupled two-leg ladders and dimerized chains in a magnetic field is discussed. These 
results are most pertinent to recent and forthcoming experiments. Finally, we examine 
weakly coupled N-leg ladders with N>2. In this case, multiple ordering transitions can 
occur which may partially overlap. In our conclusions, we propose that these overlap- 
ping high-field phases are multi-component Bose-Einstein condensates, consisting of 
magnons with different integer spin. 

2 Spin-1/2 Heisenberg Chains 

Let us consider a crystal of weakly coupled antiferromagnetic Heisenberg chains in a 
magnetic field, described by the Hamiltonian 

i/^^ = E [JiSfSf^, + SfSl, + AS^S^^,) - hS^] , (1) 

i 

where J > is an antiferromagnetic exchange constant within the chains, A is an 
easy-axis anisotropy and h is an applied external magnetic field. The chains are weakly 
coupled by H' = J' E<i,j> Si ■ Sj with < J' < J. 

The phase diagram of H^^ can be obtained from a numerical solution of the Bethe 



ansatz equations [^2|, |T3[, and is shown in Fig. 1. At zero magnetic field the ID 
subsystem is in the ferromagnetic Ising regime for A < — 1 . In the interval — 1 < 
A < 1 it is in the gapless XY regime, whereas for A > 1 it is in the massive Ising 
antiferromagnetic regime. The magnetization of the system becomes non-zero if the 
magnetic field h exceeds a minimal field given by /^mm = for — 1 < A < 1. In the 
gapped case (A = cosh 7 > 1): 
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It saturates at a maximum critical field given by /imax = (1 + A) J. 

The inter-chain coupling J' is assumed to be small compared to the intra-chain 
coupling J. Therefore, 3D long-range antiferromagnetic order only occurs in the gapless 
region. In order to calculate the magnetic response of the chains and the 3D ordering 
temperature one needs to know the low-temperature behavior of the susceptibility at 
finite magnetic fields. This can be determined analytically by mapping H^^ onto 
an effective continuum field theory that describes the low-lying excitations of H^^ in 
the gapless incommensurate regime. For A close to zero the long-wavelength limit of 
H^^ can be studied via bosonization techniques, leading to a c = 1 conformal field 
theory (CFT). The effective parameters of this CFT for the whole range of A can 




Figure 1: Phase diagram of a spin- 1/2 Heisenberg chain with an easy- axis anisotropy, 
A, in a magnetic field, h. FM: ferromagnetic regime; incommensurate: partially po- 
larized gapless regime; AFM: antiferromagnetic regime. 



be determined by comparing the thermodynamical properties of this CFT with the 
numerical thermodynamic Bethe ansatz solution to the original Hamiltonian (|I]). 

In the following, we discuss the technical aspects of this procedure. The Hamilto- 
nian (0) can be mapped onto a model of interacting spinless fermions via a Wigner- 
Jordan transformation [|l^ 
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Sf' = S'f + iSf = a} exp iir ^ a^aj , S, 



Hii H/'i 



(3) 



Using the fact that H {J, A, h) is related to H {—J, —A, h) via a unitarity transfor- 
mation, this mapping gives 
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(4) 



The magnetization of the chain, m, is related to the spinless fermion density, n , 
by m = n — 1/2. Hence, the lower critical field h^^ is determined by the condition 
that the band of spinless fermions starts to fill up. Since in this low-density limit the 
interactions between the fermions are negligible, one easily finds h^-^^ = — (1 + A) Jg. 
From a particle-hole transformation a|. — > bk, the upper critical field /i^2 = (1 + ^)Je 
(high-density limit for the spinless fermions) is obtained analogously. 

In the gapless region (i.e. hd < h < hc2) we bosonize H^'^ , and obtain a c=l CFT 
of a compactified scalar, 

Hn = I dx i 



\ 2 



-U^ + 



u 



2tiK 



{d^y 



(5) 



where (j){x, t) is the bosonic field and n(x, t) is its conjugate momentum. The Luttinger 
parameters K and u depend on the magnetic field and the exchange interaction of 
the original Hamiltonian, and still need to be determined. One can interpret Hb as 
describing a compactified boson with radius 

R=^^. (6) 

We will now discuss the equations to determine K and u, which can be derived 

from the Bethe-ansatz for H^^ (see e.g. [|15[)- One finds a system of integral equations 

1^, [l^, that need to be solved numerically. The dressed energy, ed{ri), satisfies the 



integral equation 



1 f^ 

ediv) = eo(^) - — / K{7] - r]')ediv')dv' , (7) 

Zn J-A 

where the kernel K{ri) and the bare energy eo(^) are given in Table I W7^. A similar 



table has been given by Cabra et al. [|13 
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Table I: Functions used in the integral equations for the XXZ chain. 



The cut-off parameter A is determined by the condition 

ed(A) = 0. 



(8) 



Once A has been obtained, the dressed charge function ^(r/) can be calculated from the 
integral equation 

e(^) = 1 - ;^ r i^(^ - V)av)dr^\ (9) 



(10) 



27r7-A 
directly giving the Luttinger exponent K: 

K = ^Af 



Simultaneously, the integral equations for the phase-space densities a{ri) and p{r]) are 
solved self-consistently: 
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^iv) = giv) - TT / ^(^ ~ r]')a{r]')dr]' 
Zn j-A 
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Then the spinon velocity u and the magnetization m are given by 

u e 

1 ^ 2na{Ay 
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(13) 
(14) 



where 



rfeo(A) 



^oiv)piv)dv- 



In Fig. 2 the magnetization curves, m{h), for an XY-hke, an isotropic, and an 
Ising-hke Heisenberg chain are shown, along with the upper and lower bound of the 
spinon excitation spectrum. The cases A = 0.5 and A = 1.0 are gapless, whereas for 
A = 2.0, the zero-field spin gap (Fig. 3(f)) leads to a plateau in m^h) (Fig. 3(c)). 
Furthermore, the bandwidth of the spinon spectrum increases with A, reflecting the 
renormalization of the spinon velocity due to the backscattering processes. 
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Figure 2: (a) - (c): magnetization curves, m{h), of an antiferromagnetic spin-1/2 
XXZ chain, (a) XY regime with Ising anisotropy A = 0.5, (b) isotropic Heisenberg 
point with A = 1.0, and (c) Ising regime with A = 2.0. The graphs in (d)-(f) show the 
corresponding lower and upper bounds of the spinon excitation spectrum, uj{k). 



From the exact numerical solution of the continuum model, the spin-spin correlation 
functions of the original model (|l]) can be derived. In turn, the finite-temperature 
susceptibilities are given via a Fourier transformation of the correlation functions [|T^ . 
The low-temperature behavior of the susceptibility in the gapless regime is determined 
by the dominant low-frequency spinon modes at momentum q^ = vr: 
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where B{x, y) is Euler's beta function, and F{/\) is a prefactor which strongly depends 
on the Ising anisotropy [|1^] . 
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where cos(7r/3^) = A. Here, it has been assumed that the chains are parallel to the z- 
axis of the crystal. Furthermore, we have neglected higher-order logarithmic corrections 
which arise in a more rigorous treatment of the backscattering processes. 

The low-temperature transition line to 3D ordering due to the small but finite inter- 
chain coupling J' can be calculated within an RPA approximation. The corresponding 
3D susceptibility is then given by [0, Q 

'i^(g„u; = 0;T) 



X^^(q,^ 



0;T) 



X 



(19) 



i + ^7(q)x^^(g.,^ = 0;r)' 

where /(q) is the crystal form factor which we here set to /(q) = — 1 for simplicity 
(simple cubic lattice). The 3D ordering transition is driven by the low-temperature 
divergence of x^^iQz-, ^ = 0; T), where the transition temperature is given by the locus 
of the divergence of x'^^(q, cj = 0; T). The resulting magnetic field dependence of T^ 
is shown in Fig. 3 for the various regimes of A. The field dependence of the onset of 
the crossover from ID to 3D is shown by the dashed line, which is obtained from the 
fluctuation formula derived in the next section. 

In this section, we have calculated the Luttinger parameters {K and u) of the 
effective CFT (|]) for the XXZ chain from a numerical solution of the Bethe ansatz 
equations within the whole range of A. However, the particular case of the XY limit 
(A = 0) can be treated exactly via bosonization (after linearization of the dispersion 
law around the Fermi points), since the corresponding fermionic theory is free in this 
case. In this limit, one finds the following values for u and K: 
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In this limit, the 3D transition temperature is 



T, = Jt'J/'C, 



J' VF 



K = 1. 
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(20) 



(21) 



where Ci and C2 are given numerical constants. At the lower critical field hd = — 1 
the behavior of T^ is thus T^ oc \/h — h^i |]2T|. From the numerical results a similar 
scaling behavior at the edge of the gapless phase is obtained for nonzero values of A. 
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Figure 3: 3D ordering temperature (solid line) as a function of the applied magnetic 
field in a cubic crystal of weakly coupled antiferromagnetic spin-1/2 XXZ chains, (a) 
XY regime with A = 0.5, (b) Heisenberg point (A = 1.0), and Ising regime with 
A = 2.0. The dashed lines indicate the onset of the fluctuation region below which the 
3D magnetic correlation length becomes comparable to the inter-chain spacing. For 
this plot, we have chosen J' / J = 1/16. 

3 Ginzburg Criterion 

To determine the quality of the RPA theory described in the previous section, a 
Ginzburg Criterion will now be used to examine the width AT^ of the critical re- 
gion. This can be done within the spinless fermion picture of the coupled chains, for 
which a Landau- Ginzburg functional can be derived. Denoting the localized spins by 
Sj^, where the index /x labels the ID-chains and the index i the position along these 



chains, the Hamiltonian for the entire 3D system is given by 

M i,<fM,U> 



(22) 



Introducing spinless-fermion creation operators on each chain, this Hamiltonian can 
be mapped onto a model of weakly coupled metallic chains. Neglecting the inter-chain 
hopping (assuming an Ising-like coupling between the chains), the resulting Hamilto- 



man is 
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Note that the inter-chain couphng renormahzes the bare chemical potentiaL After 
hnearization of the dispersion around the Fermi points one finds a generahzed Landau- 
Ginzburg functional, describing the 3D ordering transition of the original model in 
terms of a density-wave-type phase transition within the spinless fermion picture: 



F[^{x,y,z)] = ^ J d'x 



Ai^r + 5l^r + Ci\ V I ^r + c» ivn^i 



(24) 



Here d± denotes the interchain distance, and the other parameters can be expressed in 



terms of those of the microscopic model. Analogous to Ref. |^, they are 
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where v is the Fermi-velocity and 7 = 2 — ^^^ is the scaling exponent of the singular 
part of x^^IQz, uj = 0]T) for T —>■ 0. If T is close to T^ it follows that 



A-A'tl-. 



(29) 



with A' = J7. In the Gaussian approximation [^, the correlation length parallel to 
the chains is then given by 

eil = eon (^ - 1 j , (30) 

where the longitudinal coherence length ^o|| is defined as 






(31) 



In the direction perpendicular to the chains the correlation length is given by a similar 
expression, with 

W = /§. (32) 

Before discussing the Ginzburg criterion, let us derive the crossover condition which 
has already been mentioned above. There are different criteria for the definition of a 
crossover condition. Here the onset of the crossover is defined as the temperature, at 
which the perpendicular correlation length equals the distance between the chains, i.e. 

a = ^±- (33) 

This gives the following condition for the ID-susceptibility 

J'x^^fe,^ = 0;T) = | (34) 



The width of the critical region according to the Ginzburg criterion |^4| is given by 

AT, 1 



^> (35) 

where 



T, 32(71 AC eo±^ 



AC = [f^ (36) 

is the specific heat jump per unit volume at the transition. Using the above equations 
one finds 

^ « 0.01 1^^^^<«^4^^^. (37) 

The width of the critical region is thus of the order of one percent. 

4 Two-leg ladders and dimerized chains 

In this section, we derive effective Hamiltonians, describing the low-energy spectrum 
of various gapped spin systems that are driven into a gapless phase by an external 
magnetic field. First consider a general two-leg ladder system, described by the Hamil- 
tonian: 

H = 2_^ Jij ^i,j ' ^i+l,j ~^J±Z_^ ^i,l ■ Si,2 + Jl2 /^ ^i,l ' Si+1,2 + J2I Z_^ ^i,2 ' Sj+l,! ^ h 2_^ S^j 
i,j i i i i,j 

where 

Ja , a i even, j = 1 

J,, if i odd,, = 1 

' Jc , li i even, j = 2 
Jd , a i odd , j = 2 
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Figure 4: The two-leg ladder with various couplings. 

The various couplings are shown in Fig. 4. To derive the effective Hamiltonian in the 
gapless regime of H, a perturbation expansion in the off-rung couplings Jij and J^j is 
performed. 

Consider first the case where all couplings vanish, except for J±. If no magnetic field 
is applied, the groundstate of the ladder consists of independent singlets at each rung. 
The Hilbert-space of a single rung is spanned by the singlet | 0,0 >= (|ti> — |iT> 
)/V2, and the triplet | 1,1 >=|TT>, I 1,0 >= (|n> + |iT>)/V2, I 1,-1 >=lii>- 
The triplet excitation above the groundstate is AE = 1J± for zero magnetic field. If 
a magnetic field is applied the triplet states split, and for a critical field hcrit = AE 
the states | 0, > and | 1, 1 > become degenerate. Upon increasing the magnetic 
field further, the triplet state | 1, 1 > becomes the new groundstate of the system. 
Thus at the critical field hcrit the magnetization changes discontinuously from zero to 
saturation. 

Now consider the ladder with non-zero but small inter-rung couplings. The dis- 
continuous transition for J^ = is now broaded between the magnetic fields hd < 
hcrit < hc2, and this gapless regime is described by an effective spin-1/2 Hamiltonian. 
To derive the corresponding effective theory for this regime, one introduces effective 
spin-1/2 operators, Sf, acting on the states |lt>=| 1, 1 > and |D'>=| 0, > for each 
rung. To first order in the small couplings, one then finds an effective Hamiltonian, 
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hi" = h-J^-^{Ja + J^2 + J[2 + J2l)■ (42) 

The parameters of the effective Hamiltonian become site-independent in certain cases 
of interest. For example, if we set Jij = Jy, we recover the known result for the strongly 
coupled two-leg ladder ||25|. For Ja = Jd = J\\ = {^ ~ ^)J and Jb = Jc = 0, we obtain 



an effective description of a single chain with dimerization 5 close to one, if we set 
J± = {1 + S)J. The original model in this case is given by (see Fig. 5): 

Hd = jJ2i^ + 5i-iy)s, ■ s,,+i + j'Es. ■ s,+2 -hY^s^, (43) 

i i i 

where J' = Jjg = ^21 i^ a next nearest neighbor coupling constant. Note, that the 
isotropic chain {^eff = 1) is recovered for 



Figure 5: The dimerized chain 

In analogy to our discussion of the XXZ chain, the critical fields h^i and hc2 can 
be determined for the effective Hamiltonian, i.e. hl{ = — (1 + iS.eff)Jeff and h%^ = 
(1 + Aeff)Jeff- Expressed in terms of the physical variables of (|38D we find for the 
two-leg ladder: 

T/ I T' 

(45) 
(46) 

(47) 

(48) 

Having derived these effective Hamiltonians, the magnetic-field driven 3D ordering 
temperature can be studied for the effective models, using the results of the previous 
sections. The behavior of the transition temperature is similar to Fig. 3(c), mainly 
depending on the easy-axis anisotropy Ag// of the effective model. 
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and for the single dimerized chain: 
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5 N-leg ladders 

Let us now consider the case of A^-leg spin-1/2 Heisenberg ladders. If the couphng along 
the rungs is denoted as J^, and the coupling along the legs as J\\, the Hamiltonian is 

H^ = J|| J2 S.,r ■ S„, + J±J2 ^^,r ■ S.,.., ~hJ2Sl,, (49) 

where i and j enumerate the rungs, r, r' label the legs, and the sum marked by ^^ {\) 
runs over nearest neighbors along legs (rungs). 

We first discuss the occurance of gapless phases for these systems. Due to the 
alternating nature of the groundstates of A^-leg ladder systems at zero field, which 
show a spin gap for even N and are gapless for odd A^, one has to consider these cases 
separately. In the limit of large rung coupling (J^ ^ Jy), one can consider first a 
iV-site open Heisenberg chain in a magnetic field, and then treat the coupling along 
the leg as a perturbation. In the case of even N, there will be N/2 changes in the 
nature of the ground state upon increasing the magnetic field. At each change a new 
gapless phase opens up. This can be seen as follows: the spectrum contains multiplets 
of multiplicities m = 1, 3, 5, ..., A^ + 1. Let -^^(0) be the energy of the lowest m-plet at 
zero magnetic field. Then there is a gapless phase at a field 

A^ 

hi = E2i+l{0)-E2^-l{0), 1 = 1,...,— (50) 

In the case of odd A^, the groundstate of the system without a magnetic field is a 
doublet, and therefore the first gapless phase occurs already at hi = 0. Since the 
spectrum contains also multiplets with m = 4, 6,...,A^ + 1, there are additional gapless 

phases at fields 

A^ + 1 

hi = E2iiO)-E2^-2iO), ^ = 1,...,^— (51) 

The values of the magnetic field where the gapless phases occur can therefore be ob- 
tained by a numerical calculation of the zero-field spectrum. In each gapless phase 
there are two degenerate states, e.g. 

AT = 3, h = h2 |^>=|TTT> (52) 

W>= ^ (- ini> +2 |in> - liiT>) (53) 

For non-zero values of J\\, these two excitations spread over the ladder in the vicinity 
of the critical field, thus broadening the gapless phase. To first order in Jy, the ladder 
system can then be described by an effective spin-1/2 XXZ Heisenberg model in an 
effective magnetic field he/f, after defining effective spin operators between the states 
|f|'> and |J|> in the same matter as for the case of the two- leg ladder. Considering two 

13 



such A^-site chains, coupled to each other by a constant Jy , the energy spectrum of the 
corresponding Hamiltonian is calculated and compared to the spectrum of the effective 
model. One can thus determine the values for the parameters of the effective model 
{Jeff, ^effi heff), whcre the effective magnetic field is h^ff = h — hc{0)J± — ChJ\\- The 
obtained values are given in Tab. II. Using these effective model descriptions of the 
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Jeff/J± 


Ae// 


hM 


Ch 
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0.5 
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0.5 
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1 


1 


1 
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0.5 


1.5 


0.5 
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1 


1.0750 


0.3489 


0.6589 


0.375 
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1 


0.3750 


1.7071 


0.625 


5 


1 


1.0169 


1 
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1.0961 


0.3789 


1.1189 


0.2958 


5 


3 


1 


0.3 


1.8090 


0.7 


6 


1 


1.1114 


0.3163 


0.4916 


0.3515 


6 


2 


1.1348 


0.3011 


1.3860 


0.3985 


6 


3 


1 


0.25 


1.8660 


0.75 


7 


1 


1.0344 


1 








7 


2 


1.1415 


0.3407 


0.8848 


0.2166 


7 


3 


1.1663 


0.2381 


1.5504 


0.4891 


7 


4 


1 


0.2143 


1.9010 


0.7857 


8 


1 


1.1364 


0.3020 


0.3926 


0.3432 


8 


2 


1.1882 


0.2743 


1.1506 


0.2888 


8 


3 


1.1917 


0.1962 


1.6577 


0.5555 


8 


4 


1 


0.1875 


1.9239 


0.8125 



Table II: Parameters of the effective low-energy model for the gapless regions of N- 
leg-spin-1/2 ladders in a magnetic field, h. The effective magnetic field is given by 
heff = h- K{0)J± - ChJ\\. 

N-leg ladder in the gapless phases, we apply the RPA approach described above to a 
crystal of weakly coupled N-leg ladders. We now observe a cascade oi N/2 ((A^+l)/2 3D 
ordering transitions for quasi 1-D ladder subsystems with an even (odd) number of legs, 
as shown in Fig. 6. In the case of weakly coupled even-leg ladders, the first transition 
is driven by the formation of a SDW of triplets along the ladder direction, with a 
groundstate wavevector which is proportional to the magnetic field., h > hd- The 
following transition (for A^ > 2) is driven by a SDW of quintuplets, etc.. Depending on 
the ratio J\\/J±, these phases of different multiplet polarization may overlap, and mixed 
regimes can occur. The resulting 3D ordering temperature does not vanish completely 
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Figure 6: 3D ordering transition temperatures of N-leg spin-1/2 Heisenberg ladders as a 
function of an external magnetic field. Cascades of transitions are observed for A^ > 2, 
driven by ID SDW's of spin-S multiplets on the ladders. For this plot, we have chosen 
an anisotropy ratio J\\/ J^ = 1/5 and residual inter-ladder coupling J'/J\\ = 1/16. 

in this case, but has minima at particular magnetic fields where the number of the lower 
multiplet excitations equals the number of the next-higher multiplet excitations. The 
dependence of this overlap on the ratio of the coupling constants is shown in Fig. 7, 
where the dark stripes show the phases of 3D order at T = for the cases A^ = 3, 4, 7 
and 8. These diagrams are in good agreement with results obtained by numerical 
finite cluster diagonalizations [O]. Odd- leg ladders also have a sequence of ordering 



transitions, with the only difference that the onset of the first transition occurs already 
at h = 0. (Fig. 6 (b) and (c)). 
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Figure 7: Schematic phase diagrams of weakly coupled N-leg ladders in a magnetic 
field, h, at T = 0. The shaded areas indicate 3D ordered phases, and shown is the 
dependence on the ratio of the ladder couplings J±/J\\- 
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6 Conclusions 

We have presented a theoretical approach to magnetic-field induced 3D ordering tran- 
sitions in strongly anisotropic antiferromagnetically correlated spin- 1/2 compounds. 
These systems consist of weakly coupled chains or ladders which have singlet ground- 
states with finite excitation gaps to the lowest triplet states. Because the inter-chain 
couplings are weak compared to the intra-chain couplings, the ID subsystems are ef- 
fectively independent. By applying an external magnetic field, the spin gap in the 
subsystems can be overcome, and they become partially polarized. In this incommen- 
surate regime, even infinitesimal inter-chain coupling leads to 3D long-range ordering 
at low temperatures. 

There are only a few known realizations of Bose-Einstein condensation (BEC). The 
two most prominent examples are ultracooled dilute ensembles of trapped atoms and 
the superfiuid transition of Helium 4. In principle, BEC can occur for any bosonic 
many-body system in a confining potential at low temperatures. Therefore, magnetic 
compounds with integer-spin excitations (e.g. spin-1 magnons) present a promising 
class of candidate materials. However, indications for BEC in such materials have not 
been found until two recent experiments on TlCuCls and Cu2(C2Hi2N2)2Cl4. These 
compounds are strongly anisotropic, and have a finite spin gap. A, between their 
singlet groundstate and the first triplet excitation. Hence their magnetization is ex- 
ponentially activated at small temperatures and zero magnetic field. An applied mag- 
netic field can decrease the singlet-triplet excitation gaps of the ID subsystems, and 
eventually drive them into a partially polarized, gapless regime if the field exceeds a 
critical strength. Due to residual magnetic couphngs between the subsystems, a low- 
temperature ID to 3D transition occurs. The resulting 3D ordering may be viewed as 
a BEC of spin-1 magnons, accurately predicting the temperature dependence of the 
critical field (hBE(T) — A oc T^/^) and the magnetization curveQ, |]. 

Let us now consider multi-component phases, as they occur in spin ladders with 
more than two legs. In a spin ladder, a spin-1 magnon excitation essentially corre- 
sponds to the formation of a spin triplet on a rung, a spin-2 magnon corresponds to a 
quintuplet, etc. . At intermediate fields (larger than hd) these excitations can coex- 
ist, depending on the choice of parameters. The resulting low-temperature condensate 
may then contain multiple components, as it is the case at high fields in Figs. 6(c) and 
(d) (dashed regions). These phases may be viewed as multi-component Bose-Einstein 
condensates, consisting of magnon excitations with different integer-spins. It is thus of 
interest to conduct high magnetic field experiments on appropriate candidate materials 
to determine whether such phases exist in physical systems. 
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